Draft version September 14, 2009 

Preprint typeset using I4TgX style emulateapj v. 4/12/04 



ON THE REDSHIFT EVOLUTION OF MGII ABSOPRTION SYSTEMS 



0^ 
O 
O 
(N 

D 
00 



o 

6 

> 
On 
rn 

00 
O 
On 
O 



Jeremy L. Tinker' & Hsiao-Wen Chen^ 

'Berkeley Center for Cosmological Physics, University of California-Berkeley 
^ Department of Astronomy and Astrophysics & Kavli Institute for Cosmological Physics, University of Chicago 

Draft version September 14, 2009 

ABSTRACT 

We use a halo occupation approa ch to connect Mg II absorbers to dark matter halos as a function of red- 
shift. Using the model constructed in [Tinker & Chen (2008), we parameterize the conditional probability of an 
absorber of equivalent width Wr being produced by a halo of mass M/, at a given redshift, P(W, |M/,,z). We con- 
strain the free parameters of the model by matching the observed statistics of Mgll absorbers: the frequency 
function f{Wr), the redshift evolution n(z), and the clustering bias bw- The redshift evolution of W,- > 1 ab- 
sorbers increases from z = 0.4 to z = 2, while the total halo cross section decreases monotonically with redshift. 
This discrepancy can only be explained if the gaseous halos evolve with respect to their host halos. We make 
predictions for the clustering bias of absorbers as a function of redshift under different evolutionary scenarios, 
eg, the gas cross section per halo evolves or the halo mass scale of absorbers changes. We demonstrate that the 
relative contribution of these scenarios may be constrained by measurements of absorber clustering at z > 1 
and z ~ 0.1. If we further assume a redshift-independent mass scale for efficient shock heating of halo gas of 
^crit = 10"^ M0, absorber evolution is predominantly caused by a changing halo mass scale of absorbers. 
Our model predicts that strong absorbers always arise in ~ Man halos, independent of redshift, but the mass 
scale of weak absorbers decreases by 2 dex from < z < 2. Thus, the measured anti-correlation of clustering 
bias and Wr should flatten by z 1 .5. 

Subject headings: cosmology: observations — intergalactic medium — quasars: absorption lines — dark matter: 
halos 



1. INTRODUCTION 

Absorption lines in the spectra of quasars reflect the pres- 
ence of gas intersected along the lines of sight to the QSOs. 
A subset of these absorption lines systems, i.e., Mgll, CIV, 
Lyman-limit systems, are created when a line of sight passes 
through a dark matter halo that contains cool or warm gas. 
The presence of Mg''" indicates that the gas is photoinized with 
a temperature near lO** K. The connection between Mgll ab- 
sorbers and dark matter halos is well established: Mg II ab- 
sorbers correspo nd to neutral hydroge n column densities of 
1018-1020 cm-2 dChurchill et a l. 2000; Rao et alj|2 006), den- 
sities that exist in collapsed objects. Blind surveys of galax- 
ies near QSO sig htlines y ield Mg II absorbers a high fra ction 
of the time (Tri pp &Bowe n 2005; Chen & Tinkeiil2008l) . and 
known Mg II absorbers are nearly always found to be asso- 
ciated with a galaxy with a projected separation significantly 
less than the virial radius of the dark matter halo around that 
galax y (e.g., Bergeron 1986; Steidel 1993; Kacprzak et al. 
|2008|) . Mg II absorbers thus probe the cool gas around galax- 
ies from which these galaxies form their stars. The redshift 
evolution of absorbers is thus connected both to the growth 
of dark matter structures and the evolution of gas accretion 
in dark matter halos. Connecting quasar absorbers with dark 
matter halos is a direct probe of the gas accretion history of 
such objects. 

The observed statistics of Mgll absorbers are a conse- 
quence of the properties of the halos in which the gas resides. 
In this paper we present a model that provides a mapping be- 
tween rest-frame Mg II absorption equivalent width, Wr(2796) 
(hereafter Wr) and dark matter halo mass Mi, in such a way 
as to match the frequency function of absorbers, their large- 
scale clustering bias, and the redshift evolution of their num- 
ber counts. 

The methods for quantifying the connection of galaxies to 



dark matter halos are now well established through halo oc- 
cupat i on models (see, e.g., [S eljak^ '2000'; 'Scoccimarro et alj 
l200lt ICoorav & Shedil lIOOalB erlind & Weinberg 2002 for 
early works and Zheng et al. 2007i; ivan den Bosch et al.il2007l : 
Tinker et al. 20Q7 for examples o f more recent extensions of 
the analysis). In iTinker & ChenI (|2008) (hereafter Paper I), 
we presented a model of the halo occupation of cold gas that 
addressed the first two of these issues: the frequency and bias 
functions at z = 0.6. Our approach in Paper I was to param- 
eterize the conditional probability of an absorber being pro- 
duced by a halo of given mass P(W, |M/,) and using the statis- 
tics of dark matter halos to constrain the free parameters of 
P(Wr\Mh) in order to match the observations. After adopting 
a cored isothermal density profile for the Mg II gas (a choice 
later supported in Chen & Tinker 2008, hereafter Paper II), 
the two free quantities governing P(W, |M/,) were constrained 
from the data: the incidence of absorption per halo Kg{Mh) 
and the absorption efficiency AwiM/,). The former encapsu- 
lates both the covering fraction of Mg II gas in the dark matter 
halo and the incidence of extended gaseous halos. The latter 
controls the mean absorption strength as a function of halo 
mass. Because the halo-based model parameterizes these two 
quantities, it makes no assumptions about the physical nature 
of these systems — whether they are produced from infalling 
or outflowing gas, or within a turbulent, static gas halo. The 
primary purpose of this model is to determine the distribution 
of halo masses probed by Mg II absorbers and the incidence 
required to match observations. 

The difficulty in reconciling halo statistics with measured 
statistics of Mg II absorbers at z = 0.6 lies in the observed 
anti-correlation between W, an d the large-scale clustering bias 
(, Bouche et al. 2006: Gauthier et alJl2009t iLundgren elan 
l2009ir For dark matter halos, the most biased objects are 
the rarest and most massive. The anti-correlation between 
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W,. and bw implies that the most biased absorption systems 
are the weakest, most frequent systems. In Paper I we con- 
cluded that the fiWr) and bw data could only be fit by a 
model that contained a critical mass scale below which ab- 
sorption efficiency is high, but at higher masses the ab- 
sorption efficiency is significantly attenuated due to a lack 
of cold gas in the halo. This transition from 'cold-mode' 
gas accretion to 'hot-mode' accretion is motivated by results 
of analytic models and numerical simulations that demon- 
strate a critical halo mass at which shock heat ing becomes 
efficient (Birnboim & Dekell 1200 31; iKeres et a l.1 12005, 2008; 
Dekel & Birnboim 2006; Birnboi met all 12007; Broo ks et alJ 
2008,) . These simulations found that the transition between 
the cold mode and hot mode occurs between 10" /;~'Mq 
and lO'^ M0, precisely the location required to match the 
Mg II statistics. The expectations of these papers are qualita- 
tively similar to those from tw o-p hase gas cooling m o dels o f 
iMo & Miralda-Escudd ( Il996h and lMaller & BuUockl (|2004 . 
In these models, in halos below ~ lO" ^ H'^Mq gas cooling 
is rapid, while at higher masses the fraction of cold gas in 
the halo (in the form of pressure-supported cold clouds) de- 
creases exponentially with halo circular velocity, leaving only 
a small fraction of cold gas at M/, > lO'-' M0. These ob- 
servational trends may also be produced if the distribution of 
cold gas changes within the halo as well; the cold gas frac- 
tion may remain constant but be redistributed into a smaller 
cross section within dense gas clumps or be less extended in 
the dark matter halo. Along with a high-mass cutoff, the data 
also require a low-mass cutoff below M ~ lO" Mq due to 
the high overall bias of absorbers. The theoretical motivation 
for such a cutoff is less clear than a high-mass cutoff, and we 
will discuss both equally in this paper 

Paper I demonstrated that the anti-correlation between Wr 
and bw is produced from residual amounts of cold gas in hot- 
mode halos. This small amount of cold gas can produce only 
low equivalent width absorption in Mi, > 10'^ /i"' Mq halos, 
increasing the bias of low-W,- systems. This scenario is similar 
to the "cold streams" seen in simulations (iKeres et al.l 120051: 
iBirnboim et ani2007b . 

In Paper II we used an observational sample of galaxies in 
quasar fields to test the assumptions and results of Paper I: 
namely, that the extent of the Mg II gas radius is equal to 1 /3 

1/3 

the dark matter halo radius, thus scaling as M,/ . Brighter 
galaxies occupy more massive halos, thus the implication of 
this assumption is that the radial extent of Mg II absorption 
around galaxies scales with luminosity such that 

^^=^.*(£) • (1) 

Paper I assumed that Rg^ « 80 /z"'kpc and a « 0.42, val- 
ues obtained from the mean halo radius of galaxies de- 
rived from halo occu pation analyses of galaxy clustering at 
multiple redshifts (ie. lTinker et al.ll2007l: Ivan den Bosch et aP 
l2007tlZheng et al.ll2007h . These assumptions were in good 
agreement with the observational sample of Paper II: = 9 1 
/i"' kpc and a = 0.35. Once accounting for the intrinsic scal- 
ing between galaxy luminosity and halos radius, we found 
the incidence of absorption to be > 80% for L ~ L, galaxies, 
echoing the result from Paper I that Kg must be nearly unity at 
~ 10'^ /i"' M0 to account for the frequency of absorption at 
z-0.6. 

In this paper, we extend the results of Papers I and II to 



constrain the redshift evolution of halo occupation of the cool 
gas probed by Mg II systems. We present two complemen- 
tary models of gas-halo evolution, demonstrating that mea- 
surements of the clustering of absorbers at multiple epochs 
can provide constraints on the relative contribution of each 
evolutionary scenario. We make predictions for the evolving 
connection between absorbers and dark matter halos from a 
single model constrained to match the same critical mass scale 
between cold-mode and hot-mode halos and found in Paper I 
and predicted by the hydrodynamic simulations listed above. 

Throughout this paper we assume cosmological p arame - 
ters in agreement with the results of ISpergel et al.l (l2007h : 
n„, = 0.25, nA = 0.75, n = 0.96, h = 0.7, and erg = 0.8. All 
units assume Ho = 100 h kms"' Mpc"'. All distances are in 
comoving units unless otherwise stated. 

2. HALO VS. ABSORBER EVOLUTION 

iProchter et al.l (l200 6a*) measured the redshift evolution of 
absorber counts from Data Releas e 3 of the S l oan Dig ital 
Sky Survey. Figure [T] shows the Prochter et al. (2006a) re- 
sults broken into three bins in equivalent width. They found 
that the number of absorbers per line of sight per unit red- 
shift, n(z), increases by nearly a factor of 3 from 0.5 < z < 1 .5 
for Wr > I systems, translating to roughly a constant num- 
ber of absorbers per unit comoving distance. This is the op- 
posite to the growth of dark matter structure; as redshift in- 
creases, the amplitude of dark matter fluctuations decreases, 
reducing the number of dark matter halos. The thick gray 
curve in Figure [T^ indicates n{z) for dark matter halos with 
M> 10" /;"' M0(roughly the mass scale for L O.IL* galax- 
ies), weighted by the cross section of each halo. The halo 
number counts decrease by more than a factor of 5 over the 
same redshift range as the observations. The apparent anti- 
thetical evolution of halos and absorbers indicates that the 
gaseous content of dark matter halos must be evolving with 
cosmic time. 

A change in the Mg II gaseous cross section physically cor- 
responds to a change in the amount and/or distribution of cool, 
metal-enriched gas in the dark halo population. This could re- 
sult from a change in gas accretion onto the halo or in outflows 
from the galaxy within the halo, as well as a redistribution of 
the metal-enriched cold gas within the halo. In purely geo- 
metric terms, these changes could result in an expansion or 
contraction of the extent of the gaseous halo per object, or 
it could change the total number of halos contributing to the 
overall cross section. Because the focus of our model is to 
accurately constrain the halo distribution producing the ab- 
sorption, we parameterize the geometric consequences of an 
evolving baryonic content without making assumptions about 
the mechanisms. Thus we consider two modes of gas-halo 
evolution that can counterbalance the growth of cosmic struc- 
ture. First, the cross section per halo may be changing with 
redshift. In this evolutionary mode, the halo mass regime in 
which Mg II absorption originates is constant, but the cross 
section per halo grows at a rate commensurate with the de- 
creased halo abundance. We will refer to this scenario as "ra- 
dius evolution". In the second evolutionary mode, the extent 
of the gas halo (with respect to its dark matter halo) is constant 
with redshift, but the halo mass scale evolves. In Paper I we 
found that Kg for halos of M < 10" /i"' M0 falls off steeply, 
thus these halos contribute little to the overall cross section at 
z = 0.6. In the cold mode/hot mode paradigm, all halo gas at 
this halo mass scale should be cold, thus the decreased cross 
section is more likely to be explained by a change in the distri- 
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Fig. 1. — The redshift evolution of absorber counts. Points with eiTor bars represent the measurements of Prochter et alj <2006al) broken into tiiree bins of 
equivalent width. The dashed and dotted curves represent the radius and mass evolution models, respectively, described in §2. The solid curves represent the 
combined model, discussed in §5. The thick gray line in panel (a) shows it evolution of dark matter halo counts for M;, > 10" hr^ Mq multiplied by a factor of 
0.01 for comparison to the data. 



TABLE 1 

Parameters of the Best-Fit Models 



Model 


Aw 


CtA 




7sh 


log Men, 


logK, 


log K2 


log fl3 


logK4 


/cold 


/3m 






Mass Evol. 


23.4 


0.208 


1.57 


1.19 


11.1 


-1.60 


-0.445 


-0.053 


-0.114 


0.108 


3.51 




117.7 


Radius Evol. 


10.2 


0.047 


0.566 


0.638 


11.9 


-2.48 


-1.02 


-0.001 


-0.001 


0.057 




\A1 


122.7 


Combined 


14.5 


0.197 


1.00 


1. 00 


11.5 


-2.00 


-0.520 


-0.130 


-0.5477 


0.114 


2.52 


0.41 


122.4 



bution of cold gas or metal content than a change in the tem- 
perature structure of the gas. We discuss a possible physical 
interpretation of the behavior at low mass in §6.2. If this lower 
mass limit decreases with increasing redshift, it would com- 
pensate for the lower overall amplitude of the mass function 
at high z- We will refer to this scenario as "mass evolution". ' 

Absorber number counts, either in the form of n(z) or the 
frequency function fiW,-), cannot distinguish between these 
two scenarios. However, measurements of the large-scale bias 
factor bw can break the degeneracy and discriminate between 
radius and mass evolution. These two evolutionary processes 
are complementary but not mutually exclusive; the truth may 

' There is also the possibility that absoiption efficiency. Aw, changes with 
redshift at fixed halo mass. This is effectively similar to radius evolution; 
increasing the strength of absorption expands the radius at which a given Wr 
occurs, even though the extent of the gas radius is fixed. However, without 
any change in the gas radius relative to Shalo- this form of evolution cannot 
match the observations; the increase in the cross section at high W,- comes 
at the expense of cross section at lower W,. The total number of absorbers 
integrated over all equivalent width cannot change. Because radius evolution 
entails an increase in the cross section per halo at fixed W, , we concentrate 
on the first two scenarios. 



lie somewhere in between. In a third model, we constrain 
the relative contributions of mass and radius evolution with a 
model constrained to match the expected transition mass scale 
between cold mode and hot mode accretion found in Paper I. 
We will discuss this "combined evolution" model in detail in 
§5. 

One point to stress is that dark matter halos have comov- 
ing radii that are constant at fixed halo mass. The transition 
regions between virialization, infall, and Hubble flow are de- 
termined by the background matter density, which is constant 
in comoving units as well. In physical units, dark matter ha- 
los "shrink" as redshift increases, but their density relative to 
the background is constant. This will become important when 
discussing Kg{M,z) for low-mass halos. 

3. METHODS 

3.1. The Halo Occupation Model for Cold Gas 

To construct the gaseous halo, we follow the outline in Pa- 
per I with minor modifications that we note below. First we 
choose the density profile of the cool, clumpy gas that forms 
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part of the gaseous halo responsible for Mg II absorption. We 
assume this profile to be a cored isothermal model of the form 

Pg^s(r) oc (r'^ + aj) ' . (2) 

The absorption equivalent width is proportional to the pro- 
jected density at a given impact parameter, pb- This is 
consistent with observations that demonstrate that Mg II ab- 
sorbers seen in low-resolution spectra are broken into discrete 
components — cold clumps or clouds — when observed in high 
resolution. The total Wr is linearly proportional to the number 
of components. This leads to a relationship between impact 
parameter and Wr of 



Wripb) OC 



1 



: arctan i 



-Pi 



(3) 



where Rg is the extent of the cool-gas halo. A relationship 
of this form is in good agreement with the results of Paper II. 
The results of Paper II also indicate that the extent of the Mg II 
gaseous halo is Rg « 0.47? 200 at z = 0.6 and Uh/Rg ~ 0.2. We 
define the dark matter halo as having a mean interior density 
200 times the cosmic mean density, thus 



Ri^- 



3M,, 



1/3 



(4) 



^47rp,„A 

where p„ = Pa-it^,,, and A = 200. Recall that equation Q im- 
plies that the comoving halo radius is the same at all redshifts 
for halos of the same mass. 

We normalize the Wr(pb,Mit) relation (Eq. ID) by the equiv- 
alent width obtained by a sightline that passes through the 
center of the halo. 



WAO\Mh) 



--Aw{MH)Mli\ 



(5) 



where M12 is the halo mass in units of 10'^ M©. The func- 
tion AwiMj,) encapsulates physical properties of the gas halo, 
including the fraction of Mg II gas and the properties of the 



cold clouds or clumps within the halo. If the Mg II gas frac- 
tion were fixed, more massive halos would contain more Mg II 
gas and have a larger path length at pb = 0, thus the natural 

1/3 

scaling of W,- is to increase as M,/ . The function Aw(Mh) 
represents deviations from this natural scaling. Some disper- 
sion about W, (0|M/,) is indicated from the results of Paper II, 
but this affects Wripb) at small impact parameters where the 
cross section is too small to make a significant impact on the 
calculations. The dependence of Wr on pb and M;, is then 



WripbMh)-- 



Pl/^l+i 



: arctan i 



'4- 



pI 



Pl + al 



WriQWh) 



arctan(/?g/fl/,) 



(6) 

where the factor of 1/ arctan(/?i,/flft) is to normalize the ex- 
pression, due to the fact that equation (O does not asymptote 
to unity at pb = 0. 

Above a certain mass scale, shock heating eliminates the 
majority of cold halo gas, attenuating the absorption effi- 
ciency of the halo. We introduce a free parameter /coia to set 
the fraction of remaining cold gas once the halo becomes fully 
shock-heated. The resulting W{pb\Mh) in the hot mode is sim- 
ply equation (|6]l multiplied by /cold- In between the cold and 
hot modes, gas heating occurs as an inside-out process, re- 
sulting in a shock-heated core of gas with core radius th at in- 
creases with halo mass (see, eg. lDekel & Birnboiml2006l) . We 
parameterize the dependence of the shock radius with mass as 



^='^sh + 7shlogMi2, 



(7) 



and restrict Rsh/Rg to be no smaller then and no larger than 
1. Inside R^^, the absorption efficiency is reduced by 1 -/coW- 
The form of equation dTJi is motivated by the simulation results 
of Keres et al. (2005), in which the fraction of cold accreted 
gas decreases linearly with logM/, through the transition re- 
gion. We will refer to the "transition mass" Merit as the mass 
at which /?sh 7^^ = 0.5. 
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It is not necessary for a sightline that passes through a halo 
to produce any absorption; the clumpy nature of the cold gas 
may create a covering fraction of less than unity. Alternately, 
only a fraction of halos may contain extended Mg II gas halos. 
We define the product of the covering fraction and incidence 
as KgiMh), the probability that a halo produces absorption if 
intersected along the line of sight. We will define the func- 
tional form of Kg{Mit) in the following section. See Paper I 
for full details on the calculation of P{Wr,Mi,). 

3.2. Parameterizing the dependence ofP(Wr) on halo mass 
and redshift 

In Paper I we set Aw to be a constant with halo mass. This 
assumption was influenced more by statistics than theory; al- 
though there is no reason to expect that absorption efficiency 
should remain constant with halo mass, the data used were not 
sufficient to constrain a mass-dependent Aiy. The addition of 
niz) ameliorates this problem. We add an additional degree of 
freedom by allowing the absorption efficiency to vary in the 
cold mode: 



AwiMh) ■■ 



Ao(Mcrit 



/10'2) 



if Mh < Merit 
if M/, > Merit. 



(8) 



We will discuss the effect of this addition in more detail in 
§5.2. As in Paper I, we model the mass dependence of Kg 
non-parametrically; we set the values of Kg at four masses 
and spline interpolate between these masses: logM, = 10.0, 
11.33, 12.66, and 14.0, and refer to the values of Kg at these 
masses as k\-k^. 

To parameterize the mass evolution mode we adjust Kg as 
a function of redshift. In Paper I we found that Kg « 1 at 
M > lO'^ /i"' Mq, and it rapidly falls at lower masses. Thus, 
instead of altering the values of K1-K4, we assume that the 
shape of Kg{Mi,) is fixed at all redshifts, but the mass scale 
evolves with time. In other words. 



M,(z) = M,-,o exp[-/3M(z - zo)] , 



(9) 



where / = 1-4, and M,_o is the value at zo = 0.6 listed in the 
paragraph above. We choose an exponential form for equa- 
tion ^ because it de scribes th e growth rate o f individual ha- 
los well (IWechsler et al. 2002;lvan den Boschii2002) . Cluster- 
sized halos have high growth rates, while halos less massive 
than the Milky Way have little evolution between redshifts 1 
and 0. With equation (|9]l, a comparison can be made between 
the growth rate required to match n{z) and that of halos of 
multiple mass scales. 

In pure radius evolution, the values of M, are held constant 
at all z, but the gas radius evolves as 



Rg(MH,z) = R2oo(MH)xQA 



1+z 
1+zo 



/3« 



(10) 



Equation ( fTOb implies that the gas radius can evolve with re- 
spect to the virial radius of the dark matter A value of /3r = 1 
is equal to the gas radius being fixed in physical units, while a 
value of Pk = Q means the gas halo is fixed in comoving units, 
i.e., a fixed fraction of the dark matter halo radius. 

In both models we have ten free parameters: in both models 
there are four that govern the covering fraction (Mi -M4), two 
that control the cold-hot transition (K^^ and 7sii from equation 
|7]l, one for the cold fraction within the hot mode (/cou), and 
two that control the absorption efficiency parameter Aw(M) 
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Fig. 3. — Predictions for tlie bias of W, > 1 A absorbers for eacli of tlie 
tliree evolutionary models. The data point at z = 0.6 is the overall bias of 
absorbers from iBouche et alj 12006) . The bias curves represent the absolute 
bias with respect to the underlying dark matter distribution. 
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Fig. 4. — The mean logarithmic halo mass for three different values of Wr. 
(logM) evolves strongly with redshift for weak absorbers but is nearly flat 
for strong absorbers. At z < I, this produces the hw-W,- anti-correlation. At 
z ~ 1.3, Wr is roughly independent of redshift. At higher z, W,- and hw are 
positively correlated. 



(Ao and from equation[8]l. In each model there is one pa- 
rameter to govern redshift evolution (/?m and (Sr from equa- 
tions |9] and [TOl respectively). This is two more parameters 
than in Paper I, but we have added 57 data points for n(z) in 
three W,- ranges, for a total of 104 data points. 

3.3. Calculating Observables with the Model 

Once P{Wr\Mi,,z) is known, the redshift evolution of ab- 
sorber counts is obtained by integrating this probability dis- 
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tribution function (PDF) over the cross section- weighted halo 
mass function, i.e.. 



no 



-1/2 



(11) 



dMh 



dn 
dAh 



ag(Mh,z)P(Wr\Mh,z), 



where dn /dM/, is the halo mass function , for which we use the 
numerical results of Tinker et al.l (l2008b . The subscript z on 
dn I dMh indicates that the mass function is calculated at that 
specific redshift. The prefactor on the right hand side con- 
verts the linear density of absorbers from number per unit co- 
movi ng distance to number per unit redshift. We use the data 
from Prochter et al.l (l2006ah for n(z) in three bins in equivalent 
width: H;=[1.0,1.4] A, = [1.4, 1.8] A, and > 1.8 A. To 
estimate the errors on the data we add the Poisson fluct uations 
in quad rature with a 5% 'systematic' error. Prochter et al.l 
(l2006ah do not estimate the cosmic variance of their measure- 
ments, and it can be seen in Figure [T] that the point-to-point 
variations in n(z) can be significantly larger than the error bars 
we have estimated, suggesting that the true errors are under- 
estimated. 

The frequency function of absorbers is obtained by integrat- 
ing over both mass and redshift at a given Wi- 
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Fig. 5. — Constraints on the covering factor within Rg from the combined 
model. Inner and outer contours represent models with Ax~ < 1 and 4 with 
respect to the best fit model, respectively. The dashed curve is the best-fit 
model a.tz = 0.6, while the dotted curve represents Kg(M) at z = 2.0. At z = 2, 
the decrease in Kg above lO'- ft"' Mq, has little effect on the Mgll statistics 
because the abundance of hot-mode halos is significantly attenuated at this 
and higher redshifts. 



/(W,)= -J- J dzgiz) J dMh [j^^ ag(Mh,z)P(W,\Mh,z), 

(12) 

where cr„ = tt/?^, 



^tot = / dzgiz), 



(13) 



and g{z) is the surv ey completeness function at each red- 
shift (from IWochter et al. 2006a). As in Pa per 1, we con- 
struct the f(Wr) data from a combination of iProchter et al.l 
( I2006ah for > 1 A, and ISteidel & Sargent (EH^) for 
0.3A< Wr < l.OA, where we have adjusted the normaliza- 
tion of the " Steidel & Sargenti (119921) data to flie level of the 
IProchter et'aLi (120063 !) results at Wr = 1 A^. As with n(z), we 
estimate the errors through a quadrature addition of Poisson 
errors and a 5% systematic error bar 

The bias of absorbers is calculated in manner similar to 
f(Wr), but now each halo is weig hted by its large-sc ale bias 
b{Mk,z) (using the bias formula of [Tinker et al.ll2005l) . 



bwiz) = 



( dn 

dM 

\dM 



ag{M,z)bH(M,z)P{Wr\M,z). 



.m;Z). 

^ (14) 

When comparing to the data we set z = Zo = 0.6. iBouche et al.l 
(l2006h measure the bias of Mgll absorbers by cross- 
correlating absorbers with luminous red galaxies (LRGs) 
from the Sloan Digital Sky Survey photometric catalog. To 
compare the models to the data, we calculate the absolute bias 
of absorbers from equation (fl4l l and divide the calculation by 



^ Due to an error in the normalization of the IProchter et alj I2006af) data, 
we have increased the amplitude of /(W,) and n(z) by 30% (J. X. Prochaska, 
private communication). This adjustment brings the [Prochter et al. (2006^ 
data in to agreement with the n(z) fitting function given in .Prochter et al] 
)2006bi) . which is based on more recent results. 



the bias of LRGS at z = 0.6, bLRc = 1-81 dWhite et alj|2007l) . 
The bias data are in four bins of W^: [0.3,1.15] A, [1.15,2.0] 
A, [2.0,2.85] A, and [2.85,4.0] A. When compai'ing models to 
bias data, we take the number-weighted average bias over the 
given range in W,-- See Paper 1 for further details. 

When calculating observables, we truncate the integrals 
over mass at M/, = 10^ h~^MQ. We find that extending this 
limit to lower masses produces negligible differences in the 
model predictions, even at high redshift. We find the best- 
fit models by minimizing the total from the sum of the 
values from f(Wr), n(z), and bw- We use the Monte 
Carlo Markov Chain (MCMC; see, e.g.. lDunklev et alj|2005h 
method to minimize x^ and quantify the constraints on indi- 
vidual parameters. Results of the fitting are listed in Table 
1. 

4. RESULTS 

Figure[T]presents the best-fit models for both the mass evo- 
lution model and the radius evolution model (we will discuss 
the combined model subsequently). Both these models pro- 
duce statistically good fits, with x^/t^ = 1.30 for radius evo- 
lution and x^/'^ = 1-25 for mass evolution, with the caveat 
that the error bars are only estimates, and no covariance be- 
tween data points is accounted for. Figure |2] compares the 
model fits to the f(Wr) and bias data. The results for f(Wr) 
from each model are nearly indistinguishable, producing the 
observed exponential distribution of absorbers. In our model, 
the exponential distribution of absorber equivalent widths is 
produced primarily from the cutoff in absorber efficiency in 
the hot mode. At Wr < 1 A, there is evidence that f(Wr) 
becomes a power law ( ISteidel & Sargentlll992l: iNestor etal] 
l2005h . Above W, = 2 A, most all absorbers are produced in 
halos in the transition stage, M/, « 10"^ h~^MQ. Without a 
high-mass cutoff, /(W, ) is a power-law that overproduces the 
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abundance of high- Wr systems. Below 1 A, the abundance of 
absorbers is not attenuated by the hot mode, thus producing a 
steeper /(W, ). This is not a unique solution to the problem; 
it is straightforward to construct a mass-dependent covering 
factor function that maps the MJ^ halo mass function to the 
exp(-W,.). However, without the high-mass cutoff one cannot 
fit both the f(Wr) and clustering data at the same time. Both 
evolution scenarios produce excellent fits to bw at z = 0.6, 
yielding a strong anti-correlation with W^, as shown in Fig- 
ure 12] The amplitude of bw at zq is somewhat higher in the 
radius evolution model due to the lower value of K^^, which 
controls the transition mass Mem- The higher mass scale for 
absorbers in the radius evolution model increases the bias at 
all W,. 

As expected, both models produce nearly identical fits to 
n{z); the number counts alone cannot distinguish between 
these models. To match n{z) in the radius evolution model, 
the gas radius must inc rease relative to /?2{)0 as (1 +2)i-47±o.05 
iProchter et al.l (l2006ah find that n{z) increases roughly as 
(l+z)'-'*. Because the gas cross section scales as the 
redshift evolution of the cross section (at fixed halo mass) 
is much larger than the increase in the number counts them- 
selves. However, this is counterbalanced by the reduced fre- 
quency of dark matter halos at higher redshifts, producing an 
increase in absorber counts consistent with the data. How- 
ever, in order to match the number counts, the gas radius must 
expand in physical units to keep the total gas cross section 
constant because the number of halos is decreasing with look- 
back time. At this level of evolution, the gas radius is twice 
the virial radius of the dark matter halo at z = 2, with nearly 
unit covering factor and incidence. 

For the mass evolution model, the growth rate of the 
Kg mass scale is (3m = 3.5 ±0.2. This evolution rate is 
much larger than th at found in dark matter halos; from 
IWechsleret"an ('2002), halos with z = masses near 10 
hr^MQ have growth rates of ^ 0.6, while z = clusters have 
growth rates of ^ 1.3. This best-fit value of Pm is nearly the 
same as the evolution rate of the non-linear mass scale, which 
approximately sets the average mass scale of structures that 
are collapsing at a given epoch. The non-linear mass scale 
roughly indicates the halo mass above which the halo mass 
function experiences a Gaussian cutoff. The similarity be- 
tween the best-fit (3m and the non-linear mass scale is not sur- 
prising, given that the model is preserving the number of ab- 
sorbers in a cosmology in which the exponential cutoff in the 
mass function evolves. Adding lower mass halos to the total 
gas cross section at the same rate counterbalances the effect 
of the decreasing amplitude of fluctuations at higher z. 

One avenue for discriminating between these two models 
lies in their constraints on the transition between cold and hot 
modes. In order to fit the data, the two evolutionary models 
required the hot/cold transition scale be substantially differ- 
ent. For radius evolution, M^nt = lO"'^ h~^MQ, while mass 
evolution produces Merit = 10^' ' /!~'Mq. Hydrodynamic sim- 
ulations yield results in between these two values. 

An empirical approach for discriminating between these 
models is the evolution of their clustering. At different red- 
shifts the two evolutionary scenarios place the cold gas in 
different halos. Thus, they make distinct predictions for the 
evolution of absorber bias. In Figure |3] we show bw for ab- 
sorbers with Wr > 1 A as a functio n of redshift fo r both mod- 
els. The data point is taken from iBouche et al] (2006), and 
it has been multiplied by biRc to recover the absolute bias of 



the absorbers. In the radius evolution model bw{z) increases 
rapidly. The mean mass scale of absorbers in this model is 
nearly constant with redshift, but halos at fixed mass corre- 
spond to increasingly rare objects at higher redshift, result- 
ing in increased bias. In the mass evolution model, the mass 
scale decreases such that the bw z= 1-5 is the same as it 
is at z = 0.6. The large difference between the predictions 
of these two models demonstrates the potential of discrimi- 
nating between these two types of evolution based solely on 
their clustering, either at high redshift or at z ~ 0. Bias values 
in between these two extremes would represent a combination 
of the two evolutionary modes. The bias curves also diverge 
at low redshifts. The difference between the models at z = 0.1 
is not much larger than the z = 0.6 error bar, so bias measure- 
ments at lower redshift require more precise z = 0.6 data to 
constrain models. However, ongoing spectroscopic surveys 
should significantly enhance current z = 0.6 measurements^. 

5. PREDICTIONS OF A COMBINED MODEL 

Numerical results and analytical calculations place M„it at 
^ 2Q11.5 4 xhese results are in good agreement with the con- 
straint on Merit at z = 0.6 from Paper I. Thus, in our "com- 
bined" model mentioned in §2, we fix the hot/cold transi- 
tion to be Merit = lO"-^ /z"'M0, with a transition width of 
1 decade in mass. With these parameters fixed, we are able 
to leave both evolutionary modes as free parameters and con- 
strain their relative contribution to gas halo evolution. With 
the transition parameters fixed to those expected from simu- 
lations, both evolutionary models contribute, but mass evo- 
lution dominates dominates the overall evolution of gaseous 
halos ((3m = 2.52 ±0.24, (3r = 0.41 ±0.11). The fufl pai-ameter 
set is given in Table 1 . 

It is important to note that there are many ways to param- 
eterize this problem, and the data do not currently exist that 
afford us the freedom of making predictions without the use 
of priors on some of our parameters. Our choices are driven 
by the theoretical results discussed above. But as we demon- 
strated in §4, different choices for these parameters result in a 
different halo occupation of absorbers and result in different 
predictions for the clustering of absorbers and their connec- 
tion to galaxies of different luminosities. These are predic- 
tions that further observational results will be able to test. 

5.1. The Bias of Absorbers 

The evolution of absorber bias is shown with the solid curve 
in Figure [3] Because mass evolution dominates this model, 
the slope of bias with redshift is slightly steeper than from 
mass evolution alone. Figure |4] shows the mean logarithmic 
mass of absorbers as a function of redshift. For weak ab- 
sorbers, Wr= I A, (logM) is a strongly decreasing function 
of z. Strong absorbers, however, predominantly arise at M/, ^ 

1 /3 

Merit regardless of redshift; in the cold mode, Wr A„M,/ , 
thus the strongest absorbers arise in the most massive halos 
before the transition to shock heating. For ~ 2 A ab- 
sorbers, (logM) « 11.3 and it is nearly independent of red- 
shift atz > 0.5. 

At z = 0.6, Wr and (logM) are anti-correlated, yielding the 
bw-Wr anti-correlation measured. From Figure|4]it is apparent 

^ Baryon Oscillation Spectroscopic Survey, 

|http://co smology.lbl.gov/BOSS/ 

* W e note that the canonical value of 10'^ Mq obtained by 
IDekel & Birnboim ( 2006) specifies the mass at which the entire halo has been 
shock-heated, rather than the mass at which half the halo is within the shock 
radius, as we have defined it. 
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that the anti-correlation goes away at z = 1 .3, and at z > 2 Wr 
and (logM) (and also bw) are positively correlated. At z > 2, 
the abundance of halos above M„it is negligible, thus they no 
longer contribution to the clustering statistics of absorbers. 
The predictions for pure mass evolution are similar to those 
in Figure ID but for pure radius evolution the predictions are 
markedly different. In radius evolution, the shape of the bw- 
Wr anti-correlation is preserved at all redshifts. Thus, another 
method for discriminating between the different evolutionary 
modes is bias as a function of W,- at high redshift. 

5.2. Incidence and Covering Factor 

Figure |5] shows the constraints on the covering frac- 
tion/incidence rate of Mg II absorbers as a function of host- 
halo mass, K^. At z = 0.6, halos above 10'^ h~^MQ are 
consistent with covering factor near unity. Below this mass, 
the covering factor decreases roughly linearly with halo mass. 
The softer cutoff at low mass is somewhat different from the 
constraints on Kg{Mh) in Paper I but are consistent with the ob- 
servational results of Paper II. In Paper I, Kg foxM < 10" ^ de- 
creased rapidly, becoming essentially zero at ^ 10" M0. 
This result was driven in part by the parameterization used in 
Paper I, in which Aw was set to be a constant for halos in the 
cold mode. Here we allow Aw to vary as a power law, and 
in the best-fit model Aw ~ M^-^. Thus lower mass halos have 
somewhat lower absorptions efficiency and may have a non- 
zero covering fraction without overproducing the frequency 
of absorbers at low-W^. 

In the hot mode. Kg > 0.4 at Mi, > 10'^ /^"'Mq, the con- 
straints are poor on the covering fraction, although the l-a 
lower limit is Kg « 10%. The bias data drives the constraints 
toward higher Kg in this mass range; the models that reproduce 
the clustering anti-correlation all have Kg > 0.4. To match 
the strong clustering of W,- 1 A absorbers, as well as the 
high overall clustering of absorbers, a significant fraction of 
systems must be produced in massive halos. However, the 
current level of precision on the bias data are not sufficient 
to place tight constraints on Kg in the hot mode. A high in- 
cidence of absorption at this mass scale may seem difficult 
to reconcile with the fact that the cold fraction has been re- 
duced by ^ 90% in our model, but the model used here may 
be ove rly simplistic. In multiphase mo dels of g aseous ha- 
los (ie. lMo & Miralda-Esc ude 1996; Mai ler & Bul lock 2004), 
cold gas in hot halos extends out to the virial radius, either 
in the form of condensed clouds or cold streams. Thus, it 
is possible to match the data with a much lower Kg in the 
hot mode but now with Rg « Rim- The model predictions 
would be nearly identical for the current data, but the LRG- 
Mg II cross-correlation function would be markedly different 
at r < 1 /i"' Mpc scales. More study of the incidence and ra- 
dial extent of cold gas around clusters or LRGs is required to 
constrain the cold gas around massive objects. 

6. DISCUSSION 

We have presented a model that connects Mg II absorbers to 
dark matter halos in order to reproduce their observed statisti- 
cal properties; their clustering bias at z = 0.6, their frequency 
function, and their redshift evolution. In a model in which the 
transition between cold-gas halos and hot-gas halos is fixed 
to be 10"'^ /!~'Mq, we find that the mass scale of absorbers 
must significantly decrease with increasing redshift in order to 
match the data. Although the mass scale of strong absorbers 
varies little with redshift, the overall mass scale of absorbers 
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Fig. 6. — n{z) for the combined evolution model, extrapolated to z > 4. 
The circles are the lProchter et al.1 <2006al) data for W, > I A (z = 2.6) and 
W,- > 1.5 A (z = 4.3). The squares are the results of lJiang et alj {200% . The 
dotted line is the 'no evolution' model, in which the number of absorbers per 
comoving unit distance is constant. 



is maximal at z = 0. This upsizing is in stark contrast to the 
'downsizing' trend in the star formation of galaxies, in which 
the most massive systems form their stars early and the halo 
mass scale at which most the star formation is occurring in- 
creases with redshift. This model makes testable predictions 
of the clustering of absorbers at 1 < z < 2 and at z ^ 0.1. 

6. 1 . The evolution of absorber number counts 

The hot/cold transition is important in det ermining the 
slope o f n(z) and how that slope depends on Wr. iNestor et al.l 
(|2005|) find that the slope of n(z) depends strongly on equiv- 
alent width; for Wr ^ 0.3 A, «(z) is nearly constant, and the 
slope monotonically increases with absorption strength. The 
hot/cold transition, combined with the mass evolution sce- 
nario, offers a succinct explanation of this effect. Strong ab- 
sorbers predominantly exist in the cold mode, thus their fre- 
quency is enhanced as the fraction of absorbers contributed by 
10"-10'^ /i"' M© halos increases from z = 0.5 to 1.5. Weak 
absorbers are produced in both hot and cold halos, so the 
fraction of halos in the hot mode has less relevance. The 
frequency of 0.3 < W,- < 1 A absorbers in our combined 
mo del is nearly indepe ndent of redshift, in agreement with 
the iNestor et al.l (l2005h measurements. For pure radius evo- 
lution, the shape of n(z) is independent of W,-, lending further 
support to mass evolution as the dom inant scenario. 

Although the Proc hter et al.l ('2006a) «(z) data only extend 
to z ~ 2.2, it is interesting to extrapolate the models to higher 
redshifts. Figure |6] shows the prediction for the combined 
model at z > 2. The reduced abundance of M > 10" Mq 
dark matter halos at high redshift eventually wins out over 
either mass or radius evolution; the model peaks at z ~ 1.5, 
then turn over and decreases monotonically at higher red- 
shift. Lower-mass halos, M < 10''-'° /i-'Mo, still contain 
significant total cross section, but in our model the absorp- 
tion efficiency of these halos is too low to produce a large 
cross section at > 1 A. The results from the mass- and 
radius-evolution models are quantitatively similar. The dotted 
curve in Figure |6] is commonly referred to as a model with 
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redshift redshift 

Fig. 7. — Panel (a): Redshift evolution of the total gas cross section, Ktot = f^g (Rg/Rioo)^^ for the best-fit combined model. The black curves represent different 
halo masses. The key indicated logM for each line type. The two gray curves are the toy model results, discussed in §6.2, for logM = 1 1 .0 and 10.5 (upper and 
lower curves, respectively). Panel (b): Same as panel a, but results are shown for the 2a upper limit on ki shown in Figure[5] Gray curves have been shifted up 
to compare with the halo model results. 



'no evolution', meaning that the comoving number density of 
Mg Il-absorbing objects is constant and that the gas cross sec- 
tion per object is constant in physical units. For this model, 
n(z)oc(l+z)2[17„(l-i-z)3-i-17A]-i/2.|JimgitaD(|2003)have an 
estimate of n(z) at z = 2.6 (W,- > 1) and z = 4.3 (Wr > 1 .5) from 
a sample of six z 6 quasars with a total redshift pathlength of 
Az = 0.4. These data are preliminary; the mean number den- 
sity of Mg II absorbers may change when surveys of a larger 
sample of QSO sightlines become available. The current data 
are consistent with the no-evolution model. But from Fig- 
ure |6] it is clear that the no-evolution model corresponds to 
significant evolution with respect to the underlying halo pop- 
ulation, both in terms of the typical halo mass and the size 
of the gas halo with respect to the dark matter halo. When 
extrapolating the halo model, the predictions are an order of 
magnitude below the measurements at z = 4.3. Removing the 
critical mass scale and allowing all halos to have high absorp- 
tion efficiency at this redshift (as found in the analytic model 
of Dekel & Birnboim 2006 and the hydrodynamic simulations 
of lKere s et al. 2008) will not alleviate this discrepancy due to 
the paucity of halos above M^it at this epoch. 

There are two major possibilities for resolving this discrep- 
ancy; the Mg II cross-section around galaxy-type dark mat- 
ter halos increases significantly beyond their virial radius, or 
low-mass (sometimes called 'minihalos') have enhanced ab- 
sorption efficiency at high redshift. The former scenario could 
be the result of star formation in typical high-redshift galax- 
ies, where metal-enriched gas is efficiently dispersed into the 
IGM to increase the cross-sect ion of absorbing ga s. The lat- 
ter scenario is the proposal of I Abel & Mol (Il998h to explain 
the high frequency of Lyman-limit systems at z > 2. Star for- 
mation in these halos is not possible after reionization, but 
before that time it may be possible to form small amounts of 
stars and enrich the minihalos. Further data on the equivalent 
width distribution and absorber-galaxy pairs will be able to 
discriminate between these two models. 



6.2. Evolution of absorbers in the cold mode 

Our combined model predicts that the overall halo mass 
scale of Mgll absorbers reduces as redshift increases. This is 
caused both by the loss of massive halos as redshift increases 
and the the increased gaseous cross section of low-mass halos 
in the cold mode. The rate of change of the halo mass scale for 
low-W, systems is significantly faster than the growth rate of 
low-mass halos themselves; i.e., a 10" Mq halo at z = 1.5 
that produces an absorber is much less likely to produce an 
absorber when that halo evolves to z = even though that halo 
is, on average, twice as massive. 

Figure]?^ shows the evolution in the mean covering fraction 
within halos relative to the total halo cross section, Ktot = ^gX 
(Rg/Rioo)^, for the combined model. At 10'^ Ii'^Mq, this 
quantity does not change much with lookback time, but at 
lower masses the cross section at z = is below 1 % and the 
evolution with redshift is steep. 

We can try to understand the mass and redshift trends seen 
in Figure |7] through a simple picture of gas accretion. The 
evolution of Kg occurs in the cold mode — halos less massive 
than Ml, lO"'^ h~^MQ. At all mass scales, the accretion 
of baryonic matter is expected to increase with redshift. A 
convenient fitting formula is 

M^«6.6MiV^(l+z)^'^^ (15) 

(iDekel et al.ll2008h . Since there is no hot gas halo at these 
masses, and all gas is accreted cold, the infalling gas will sim- 
ply fall to the center of the halo. We assume that it takes a 
free-fall time, tff, for the gas to reach the galaxy at the center 
of the halo — although the gas is accreted with non-zero ve- 
locity, it will also have angular momentum that must be shed 
before it can reach the galaxy. The important assumption is 
that the infall time is related to the dynamical time of the halo. 
The total amount of gas in a cold-mode halo is Mg = t/j-Mg. 
Using 
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^lA(l+zT^'^Gyx, (16) 

the total halo gas scales as Mg ^ (1 +z)"^^M,|^^. The free 
fall time does not depend on halo mass because all halos have 
the same mean density of 200 times the background, which 
scales as {l+zy". We are less interested in the absolute scaling 
of the covering factor in this toy model than the scaling with 
mass and redshift, so we can normalize this relation to give 
the proper covering fraction of 10" /i"' M© halos at z = 0.5, 
Ktot = 0.02. Provided that the mean properties of the clumps — 
their sizes and densities — stay relatively constant, the cover- 
ing factor of the accreted gas should scale with the amount 
of material and the volume throughout which they are dis- 
tributed. 

Not only does the gas accretion rate increase with red- 
shift, but the physical volume in which they are spread is 
decreasing — ^recall that halos at fixed mass have constant co- 
moving radii. Thus the covering factor scales as 

Ktot ^ MgiMh , z)^20o(M/, , 

~Mf(l+z)2", (17) 

where 7?200 is now specified in physical units. This scaling 
has the proper behavior required to match the absorber statis- 
tics considered here: in the cold mode, the covering factor 
decreases with decreasing mass and increases with redshift. 
The two gray lines in Figure |7] show this scaling relation for 
lO" and 10'"^ /i~'Mq halos, once again normalized to the 
model value for 10" Mq halos at z = 0.5. For the best-fit 
model, the mass and redshift dependence is stronger than in 
our toy model. 

The variables of the model that control the evolution of Ktot 
in the cold mode are k\, K2, which set the slope of 
function of halo mass, and (3m- Due to the low contribution 
to the overall gas cross section from these masses, we are not 
able to constrain ki and K2 tightly. Figure]?}? shows Ktot(z) for 
the same halo masses but for the 2a upper limit on ki shown 
in Figure |5] In this model. Kg scales as M°-'' at low masses 
rather than M/, as in the best-fit model. In this realization of 
our halo model, the scalings of Ktot are in good agreement 
with the expectations from the toy model. To better constrain 
these parameters, tighter bias measurements are required — or 
bias measurements at higher redshift — as well as more precise 
measurements of fiWr) at Wr < 1 A. 

The presence of an absorber in M/, < Ma-it halos is indeli- 
bly linked to the rate of gas accretion. Even if absorbers are 
caused by active star formation episodes, such episodes can- 
not occur if there is no fuel coming into the galaxy. If Mg 11 
absorbers are related to star-forming objects it is tempting link 
the decrease in the cosmic star formation rate from z = 1 to 
z = with the measured n(z) for absorbers. Certainly this is the 
case to some extent; the redshift evolution of Mg 11 systems 
is strongest for the highest equivalent widths. But the lack 
of significant evolution in W,- < 1 A absorbers argues against 
such a link for all systems. In addition, n(z) is a cross-section 
weighted statistic and some care is required for quantitative 
comparison with the observed star formation rate density. 

6.3. Cold gas in hot halos 



In the combined model, fcou is ^ 10%. Nominally, this 
parameter relates to the fraction of cold gas that remains af- 
ter the halo has been shock heated. In addition to the bias 
anti-correlation, there are several additional pieces of ob- 
servational evidence that point to cold gas existing in ex- 
tended hot gas halos. First, the incidence of Mg 11 absorbers 
along lines of sight near z ~ 0.7 g alaxy clusters is s ignifi- 
cantly enhanced re lative to the field (iLopez et al.ll2008l ). The 
iLopez et"an (l2008h results favor enhancement of stronger ab- 
sorbers (Wr ~ 2) over weaker absorbers, but absorption of 
any kind reveals the presence of cool gas within the clus- 
ter virial radius. Second, the r ecent measurements o f the 
LRG-Mg 11 cr oss-correlation by ^Lund gren et al.l (12009 ) and 
iGauthier et al.l (|2009) extend the previous measurement to 
small scales — r < I Mpc. Pairs at this scale can only 
be created if LRGs and absorbers exist within the same ha- 
los; ie, the "one-halo" term from halo models of cluster- 
ing. Halo occupation analysis of the LRG autocorrelation 
function constrains LRG halo masses to be_> lO'^ h'^Mp, 
(iPadmanabhan et all 120081: iBlake et al.1 2008; Zhena et alJ 
|2008|). The multiphase cooling model of Mailer & BullocJ 
(|2P04) predicts that cold clouds are still able to form and sur- 
vive at halos with virial velocities of 300-400 kms"', with a 
cold fraction of 10% at ~ 320 km s~' . High-reso lution hydro- 
dynamic simulations by Kauf mann et alj (l2008h support the 
scenario in which cloud formation is efficient in hot-gas halos 
with cored density profiles, consistent with t he density pro- 
file of hot-mode halos in the lKeres et al.l (l2008i) results, which 
have lower resolution but cosmological volumes. 

Detailed analysis of the cross-correlation function of LRGs 
and Mg 11 absorbers at both large and small scales can con- 
strain both the covering fraction and density profile of cold 
gas in hot-mode halos, and this is the subject of ongoing study 
(J.-R. Gauthieret. al., in preparation). 

7. SUMMARY 

In this paper, we have used the frequency and abundance 
of dark matter halos to construct a model for the frequency, 
abundance, and clustering of Mgll absorption systems. One 
of the main precepts on which our model is based is the 
idea of a critical mass scale, below which absorption is high- 
efficiency and above which absorption is low-efficiency but 
non-zero. The free parameters of this model are constrained to 
match the known absorber statistics, allowing us to constrain 
the relationship between absorption equivalent width and halo 
mass. Our main findings are: 

(1) In order to match the increasing frequency of absorbers 
with redshift, n(z), the connection between absorbers and dark 
matter halos must evolve significantly from z = 0.5 to z = 2. 
Either gaseous halos must expand with respect to their dark 
matter radii, or the halo mass range probed by absorbers must 
expand to account for the data. Different evolutionary scenar- 
ios make distinct predictions for the clustering of absorbers at 
both z = and z > 1 . 

(2) If Merit remains constant at lO"'^ /;"' Mq at all redshifts, 
the dominant evolutionary mechanism is in the mass of a typ- 
ical absorber. Strong absorbers are preferentially created in 
Merit halos at all redshifts, while weak absorbers (W,- ^ 1 A) 
evolve from occupying 10'^ h~^MQ halos at z < 0.5 to 10" 
/i"'M0 halos atz> 1.5. 

(3) In our best-fit model, the cross section of M < 10" in- 
creases significantly from z = to z = 2. The scaling of Kg(M) 
with mass and redshift is consistent with a simple picture in 
which all gas is accreted cold and sinks to the center of the 
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halo in a freefall time. The increased covering fraction results 
from higher gas accretion rates and smaller volumes of halos 
at high redshift. 

(4) At z = 4, our fiducial model is not able to reproduce the 
preliminary measurements of the high frequency of absorbers. 
If these data do not change with greater statistics, this implies 
that the absorption efficiency of low-mass M < 10'" /!-' Mq 
is enhanced at these redshifts, or the cross-section around typ- 
ical galactic halos is expanded beyond their virial radii. 
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